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1 Introduction 

As the LHC started colliding protons with unprecedented beam energy, interest has risen 
in topics that had not received much attention previously, with the aim of uncovering new 
physics signals. Of these topics, that have recently seen substantial development, is the 
substructure of “fat” jets originating from the almost-collinear decay products of heavy 
resonances that are highly boosted (see for example refs. [1-13]). Many substructure 
techniques, such as filtering [1], pruning [4] and trimming [14], have been developed for 
the purpose of improving the discrimination of signals from QCD background. The lat¬ 
ter substructure techniques aid in providing cleaner and more accurate measurements of 
the properties of these resonances through: first, identifying the origin of the jet (decayed 
massive particle — signal — or plain QCD radiation — background—). Second, mitigat¬ 
ing away the jet constituent particles that have most likely originated from initial-state 
radiation, underlying-event and pile-up. 
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These techniques require in many situations calculations of QCD observables (e.g., jet 
mass, jet shapes, etc) which need special attention particularly in the vicinity of the thresh¬ 
old limit where they become highly contaminated with perturbative large logarithms as well 
as non-perturbative corrections. There has recently been some analytical work on a handful 
of the said substructure techniques with the aim to pave the way to a better understand¬ 
ing of their analytical properties (see [15] and references therein). Nonetheless and for the 
majority of QCD observables and substructure techniques the only other option available 
is resorting to numerical simulations which are based on Monte Carlo (MC) integration 
methods, and which use several approximations, e.g., Herwig [16, 17], Pythia [18, 19] and 
Sherpa [20]. These MC event generators have been very successful in describing collider 
data and are commonly used in the extraction of crucial information to boost the search 
for new physics. 

An important issue that needs addressing is the accuracy of the said MC algorithms 
and the range of validity of the approximations used therein. For instance, amongst the 
widely adopted approximations in the said MC generators is that of large- N c limit (with 
N c the dimension of SU(1V C ) group). The latter limit, which corresponds to neglecting non- 
planar Feynman diagrams, greatly simplifies the otherwise tremendously complex colour 
structure, especially at high multiplicities. However, MC generators are generally tuned 
with data from collider experiments for parameters that account for non-perturbative effects 
such as hadronisation, underlying event, etc. The process of tuning is itself vulnerable 
to erroneously ascribing neglected perturbative (observable-dependent) components, which 
might be originating from finite- N c corrections, to universal non-perturbative parameters. 
This could then potentially be a source of major discrepancy between the data and the 
predictions by the MC generators. It is thus of great importance to assess the validity 
of these approximations and make sure that neglected terms would not affect precision 
measurements. 

Amongst the issues that MC generators are meant to tackle is that of the resummation 
of large logarithms typically inherent in the distributions of most observables. These large 
logarithms are a manifestation of the miscancellation of infrared/collinear singularities at 
the matrix-element level, due to the exclusion of real-emission events in certain regions of 
phase space. For several observables of sufficiently inclusive nature, 1 i.e., global observables, 
the resummation of these logarithms is relatively straightforward and has even been achieved 
analytically to NNNLL accuracy [21]. In fact semi-numerical programs have been developed 
with the power of resumming a wide range of global observables up to NLL (CAESAR [22]) and 
even to NNLL recently (ARES [23]). However the extension of the resummation programme, 
up NNLL or even to just NLL accuracy in some cases, has seen slow progress for another 
class of observables, namely non-global observables [24, 25]. 

Non-global observables are observables that are sensitive to emissions in restricted 
angular regions of the phase space. The distributions of such observables contain logarithms 
(named non-global logarithms (NGLs)) of the scales present in the process. For instance, 

1 By “sufficiently inclusive” one means observables that are inclusive over emissions in the entire angular 
phase space. 
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the hemisphere mass distribution contains logarithms of the ratio Q/(Qp), where Q is 
the hard scale of the process and p is the normalised hemisphere mass squared. In the 
region where p <C 1, these NGLs can form large contributions to the said distributions, and 
should thus be resummed to all-orders. Up to very recently, their resummation was only 
possible numerically at large N c by means of: an MC program [24, 25] or solutions to the 
non-linear integro-differential Banfi-Marchesini-Smye (BMS) evolution equation [26]. The 
large-A^ c approximation significantly simplifies the colour flow in multiple gluon branchings 
enabling the possibility of the resummation of NGLs at least numerically. Much effort 
has recently been advocated to achieving numerical (analytical) resummation of NGLs at 
finite (large) N c . The work of Hatta and Ueda [27] exploits the suggestion of Weigert [28] 
to use an analogy between the resummation of small-Bjorken-x (BFKL) logarithms and 
that of NGLs at finite N c in a numerical fashion. They have noticed that the neglected 
finite-7V C corrections are indeed negligible in the context of e + e _ —>■ di-jets. They have 
however speculated that the situation may be drastically different for hadronic collisions. 
Furthermore, Rubin [29] numerically computed the NGLs series for both filtered Higgs-jet 
mass as well as interjet energy flow observables up to five- and six-loops, respectively, at 
large N c . In the same limit, Schwartz and Zhu [30] worked on the analytical solution to the 
BMS equation by means of an iterative series-solution up to five-loops. 

The major hindrance that one inevitably faces when attempting to compute NGLs an¬ 
alytically at finite N c is twofold. Firstly, the colour topology of a multi-gluon event requires 
evaluations of non-trivial traces of colour matrices in SU(IV C ), which become increasingly 
cumbersome starting from four-loops. Secondly and not less important, the non-Abelian 
gluon branchings increase the number of Feynman diagrams factorially at each escalating 
order to the extent that an automated way of accounting for all possible branchings be¬ 
comes inescapable. 2 Besides, there is also the issue of the various possible real, virtual and 
real-virtual gluon configurations that are eventually responsible for the miscancellation of 
soft singularities, thus leading to the appearance of large logarithms. These difficulties may 
have been the main reason for the slow progress in the resummation of NGLs at finite N c . 

In this paper we overcome the above-mentioned difficulties and present the first an¬ 
alytical calculation of NGLs at finite N c beyond leading order. Working in the eikonal 
approximation [31—36] for soft (strongly) energy-ordered partons, the first problem, i.e, 
that of colour structure, is resolved via the use of the Mathematica package ColorMath 
developed by Sjodahl [37, 38]. The latter program performs the summation of SU(IV C ) 
colour matrices in an automated way at any loop order. For the second obstacle we de¬ 
veloped a Mathematica code that accounts for all possible gluon branchings (and thus for 
all possible antenna functions) in an automated way. 3 Consequently we have been able to 
analytically calculate all squared amplitudes for the emission of soft energy-ordered gluons 
(for all possible real, virtual and real-virtual configurations) in the eikonal approximation 
fully (at finite N c ) up to five-loops. We leave the presentation of these squared-amplitudes 

2 The number of cut diagrams to consider at n loops is formally ((n+ 1)!) 2 for real gluon emission. This 
number is slightly reduced by considering on-shell particles and exploiting available symmetries. 

3 This code will be improved, in the near future [39], into a full Mathematica package capable of analyt¬ 
ically computing QCD eikonal amplitudes at (theoretically) any loop order. 
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and the method of calculation to a forthcoming paper [39]. 

With these squared-amplitudes at hand, we provide in this paper a calculation of 
NGLs at finite N c to single logarithmic accuracy for single-hemisphere mass distribution in 
e + e - —>• di-jets up to five-loops. While our calculation is full at four-loops it is incomplete 
at five-loops due to missing terms for which the squared amplitudes are not so simple to 
simplify and/or integrate. We find that the aforementioned distribution exhibits a pattern 
of exponentiation both for global and non-global logarithms. We consequently write the 
all-orders distribution as a product of two exponentials; the first being the usual Sudakov 
form factor and the second represents the “resummed form factor” for NGLs. For the sake 
of cross-checking we take the large-iV c limit of our result and compare it with previous 
calculations obtained by Schwartz and Zhu [30]. We find complete agreement up to our 
accuracy, which is four-loops. Furthermore, we compare our analytical resummed factor to 
the available all-orders numerical results [24] and discuss the phenomenological implications 
of our findings, particularly the issue of the accuracy of the large- N c limit, by assessing the 
importance of neglected finite-Ay corrections up to four-loops. 

This paper is organised as follows. In the next section we outline the usual procedure 
of calculating NGLs by defining the observable, kinematics and the general relation for 
the hemisphere mass distribution in terms of the squared amplitudes and a “measurement 
operator”. 4 We present, in the same section, the calculation of NGLs at leading order 
(two-loops) to warm up for higher loops. In section 3 we explicitly calculate NGLs beyond 
leading order at three, four and five-loops. The difficulties associated with calculations 
at five-loops will be addressed therein. We compare our findings, in section 4, to those 
obtained at large N c in ref. [30] as well as to the all-orders parametrised form reported by 
Dasgupta and Salam, which they obtained by fitting to the output of their MC program 
[24]. We also assess the relative size of the corrections due to finite N c up to four-loops and 
discuss our findings in the same section. Finally, we conclude our work in section 5. 

2 Hemisphere mass distribution at one and two-loops 

Our aim in this paper is to calculate NGLs at finite N c to single logarithmic accuracy up to 
the fifth order in the strong coupling a s (or equivalently up to five-loops). Our calculation 
is performed using QCD squared-amplitudes for the emission of energy-ordered gluons in 
the eikonal approximation. The latter is sufficient to capture all single logarithms a”L n , 
with L being the large NGL. As stated in the introduction, we do not show herein explicit 
formulae for the said squared-amplitudes and refer the reader to our coming paper [39]. 
Moreover, for the purpose of this paper, we do not consider the role of any jet algorithm, 
and postpone such work to future publications. 

2.1 Observable definition and kinematics 

For the sake of illustration and to avoid unnecessary complications from a hadronic environ¬ 
ment, we choose to work with the same observable that was used in the original paper on 

4 The idea of the “measurement operator” was introduced by Schwartz and Zhu in their paper [30] which 
we found very helpful in organizing the real-virtual contributions to NGLs. 
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NGLs by Dasgupta and Salam [24] within the framework of QCD, that is, the hemisphere 
mass distribution in e + e^ —>• di-jets. This very observable was also considered in refs. 
[40-42] in the context of soft-collinear effective theory (SCET). In both refs. [24, 40-42] 
NGLs were only computed up to two-loops. In the eikonal approximation, sufficient for our 
purpose, we consider energy-ordered soft gluon emissions: 5 Q kti 2> kt 2 2> ■ ■ • kt n , 
with Q the centre of mass energy and ku the transverse momenta of emitted gluons ki. We 
note that gluon decay into quarks has a sub-leading contribution to NGLs as was found at 
two-loops in refs. [40, 43]. 


The four-momenta of the outgoing quark, anti-quark and gluons 
pidity parametrisation as: 

are expressed in ra- 

Pq = |(1,0,0,1), 

(2.1a) 

Pq= f (1,0, 0,-1), 

(2.1b) 

ki =ku (cosh rp , cos fa , sin fa , sinh ip ), 

(2.1c) 


where recoil effects are negligible to single logarithmic accuracy. Here pt and fa are the 
rapidity and azimuthal angle of the z th emission and ku its transverse momentum with 
respect to the z-axis, which we choose to coincide with the outgoing quark direction. We 
have ku = Wj sin 9i , with cj* the energy of gluon ki and 6i its polar angle. The rapidity is 
related to the polar angle 9i through the relation pi = — lntan(0j/2). 

We define the following “antenna” functions relevant to the squared amplitudes that we 
use in this paper: 

i ,2 Pa-Pb 

Wab U {Pa-ki) (ki.p b ) 

A ab = Kb {Ki + Kb - Kb) > ( 2 - 2b ) 

^b=Kb{K!: + Kb-K^- (2.2c) 

The quark and anti-quark directions define two coaxial back-to-back hemispheres (Hl and 
Hr) whose axis coincides with the thrust axis at single logarithmic accuracy (see figure 1). 
We pick for measurement the hemisphere pointing in the positive z-axis (quark direction). 
The normalised hemisphere mass (squared) p is then defined by: 


P — | Pq T ^ ] ki J /Q 
i£H R ) 

Pi= 2 ki- Pq/Q 2 = Xi e~ Vi 


2 'y ^ ki-Pq/Q — > ] pii 

i&'H.R i 


(2.3) 


where we introduced the transverse momenta fractions Xi = ku / Q , and the sum over the 
index i extends over all emitted real gluons in the measured hemisphere Hr. 

5 Since gluons must satisfy Bose statistics one should normally allow for the permutations of the gluons 
and divide by a factor n\. This is however equivalent to choosing a specific ordering and removing the 1/n! 
factor. 
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Hr 

Measured Hemisphere 


Figure 1. Schematic diagram for an outgoing qq pair associated with multiple gluon emission. The 
measured hemisphere is the one pointing in the quark direction (Hr). 


We compute the integrated hemisphere mass distribution (cross-section) normalised to 
the Born cross-section, defined by: 



p l da 
do dp' 


— 1 + £l(p) + X 2 {P) + ' ' ‘ ! 


(2.4) 


with 


>X2>--->Xm 

where W\2—m = W x (&q, k> 2 , ■ ■ ■ , k m ) is the eikonal matrix-element squared for the emission 
of m energy-ordered soft gluons of configuration X off the primary qq pair at m th order, 
normalised to the Born cross-section. The sum over X extends over all possible real (R) 
and/or virtual (V) configurations of all the gluons {kj}. For instance, at 2 loops (m = 2) 
the eikonal squared-amplitudes Wj/, over which the sum is taken are: W R 2 R , W R r/, W^R, and 
W 12 ■ The quantity Wf/, for example, is read as: the squared amplitude for the emission of 
two energy-ordered gluons, k± and k 2 , with gluon k± real and gluon k 2 virtual. Notice that, 
in the eikonal approximation, the squared amplitude for the softest gluon being virtual is 
simply minus the squared amplitude for it being real. In other words: 

, ( 2 - 6 ) 


SmW =E / 


n 

i =1 


d$i 


U rn W X 


12 ■■■m i 


(2.5) 


where x could either be R or V. At one- and two-loops, for example, one has: 

Wi = — VVf , 

Wfv = _ yyRR ; W/v = _^vn 

The phase space factor for the emission of m gluons is: 


n«<=n 


d 3 kj 

(27t) 3 2c jj 



(2.7) 


( 2 . 8 ) 
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where g s = \fATra s and a s = ol s /tt. The factor n"=i ^ti/^ds multiplies the squared ampli¬ 
tude >Vf 2 ... m to produce the purely angular squared-amplitude (i.e., depends 

only on ?/ and 0 variables and the corresponding colour factor). In other words, one may 
write: 


where 


n d ^wf 2 ... m =n d ^wf 2 ... m 

2=1 2=1 


m m ~ o m , , , 


2=1 


2=1 

= >vf 


2=1 


TT — 

iJ -2. 9 2 
4 = 1 


(2.9) 


( 2 . 10 ) 


The non-linear “measurement operator” U m acts on the squared amplitudes and 

plays the role of excluding gluon emission events for which the hemisphere mass is greater 
than p. It is not, however, equivalent to a simple heaviside step function 0(p — Yli x i e ~ Vi )^ 
since it requires non-numerical input (information about the real-virtual nature of the 
various gluons). Due to strong ordering, the measurement operator U m factorises into a 
product of individual measurement operators; U m = ri-=i^*- The squared amplitudes 
Wi 2 —m are eigenfunctions of the measurement operators Hi with eigenvalues 0 or 1 such 
that: 


• if gluon ki is virtual then Ui ... TO > 

• if gluon ki is real and outside 'Hr then rq ... m = Wj K 2 ... m , 

• if gluon ki is real and inside Hr with Pi < p then Ui VV^... m = , 

• if gluon ki is real and inside Hr with Pi> p then Ui = 0. 

This means that events with real emissions inside the measured hemisphere and which 
contribute more than p to the hemisphere mass are excluded (i.e., not integrated over). 
That is, E(p) represents the probability that the measured hemisphere mass be less than 
p, as is expressed in eq. (2.4). We therefore write the measurement operator as: 

in =QJ + 0f [0° ut + 0“© (P - Pi)] = 1 - 0f0j n 0f, (2.11) 

with 0^ = @(pi — p) = Q(xie~ Tli — p ), 0) n = @(rn ), and 0° ut = @(—rji ). The heaviside 

step functions 0) n and 0° ut respectively indicate whether gluon ki is inside or outside the 

measured hemisphere region. The operator 0? (07) equals 1 if gluon ki is real (virtual) 
and 0 otherwise. If gluon ki is real then 0fW^... m = and 07Wj \.. m = 0, and vice 

versa. In the above we used the relations ©) n + 0° ut = 1 and Qf + 07 = 1. 
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2.2 One-loop calculation and the Sudakov exponentiation 


Having properly set up the form of the integrated distribution, we can now proceed with 
the calculation at each loop level. To warm up for higher loops we start with the one-loop 
case. At one-loop the squared amplitude for the emission of a single real gluon, multiplied 
by the corresponding phase space, is given by: 

d<f>i x W? =d s — dr]i^p- x C F w l (M 

=a s —±d m ^2C F . ( 2 . 12 ) 

X\ ZTT 

The corresponding virtual contribution is = —Wf. The measurement operator reads: 

tii = ui = 1 - , (2.13) 


which when acting on the squared eikonal amplitudes yields: 

u\ Wj 1 + uiW? = Wf - ©?0i n 0fWf + >% v - 0?0i n 0fWr 

= -0( , 0 1 1 n W 1 R , (2.14) 

where we used © R W^ = 0, 0 R Wf = Wf and Wf + Wi = 0 (which means that the 
real and virtual contributions completely cancel out in sufficiently inclusive cross-sections). 
Substituting into the expression of (e.q. (2.5) with m = 1), we are left with the 

uncancelled integration: 


£i(p) 


J d$i©?©i n Wf 


—2 CfUs 



dxi f 27T d(f>i 

xi J 0 27r 


(2.15) 


where the step function 0(* restricts x\ to be greater than p e m , and since x\ < 1 then one 
also has the restriction on the rapidity such that rj\ < L, with L = ln(l/p). Performing the 
integration to single logarithmic accuracy we find: 


S F (p) = -C F a s L 2 = ^(p). 


(2.16) 


We note that the leading logarithms in the hemisphere mass distribution are double log¬ 
arithms, which originate from soft and collinear (to the direction of the outgoing quark) 
singularities of the squared amplitudes for primary gluon emissions off the initiating hard qq 
pair. It has long been known that the resummed distribution accounting for these primary 
emissions (or global logarithms) to all-orders is entirely generated from the leading-order 
result by simple exponentiation (Sudakov form factor). However, Dasgupta and Salam [24] 
showed that a new class of large single logarithms appears starting at two gluons emis¬ 
sion, and which they termed NGLs. We may thus express the resummed hemisphere mass 
distribution as follows: 


£(p) = £ p (p)x£ NG (p) 


(2.17) 


where £ p is the primary Sudakov form factor, 

£ p (p) = 1 + £ p + i (£ p ) 2 + i (S p ) 3 + ■ ■ ■ 

= exp (X p ) = exp (—Cfc^A 2 ) , (2.18) 

and £ NG is the resummed non-global factor, 

£ NG (p) = 1 + S* G (p) + ^ G (p) + ■■■ . (2.19) 

Our aim in this paper is to compute the non-global functions X^ G (p) for m = 2, 3,4 and 5. 

2.3 Two-loops calculation and non-global logarithms 

The various squared-amplitudes for the emission of two energy-ordered (real or virtual) 


gluons with x 2 <C x\ <C 1 are expressed as (see for instance [44]): 

Wf 2 R =Wf W 2 R + Wi 2 R , Wf 2 v = - Wf 2 R , (2.20a) 

W^2 = ~ WfW 2 R , WX 2 = ->% v 2 r , (2.20b) 

where the irreducible term W R 2 is: 

W™ = \c F C a A%, (2.21) 

and the measurement operator in this case is given by: 

Uo = uxu 2 = (1 - ©?©i n ©?) (1 - 0 2 © 2 n @ 2 ) 

=1 - ©?©“©? - 0^0 2 0 2 n 0 2 (©^ + 0? ut ©f) , (2.22) 

where 0 2 = 0 P 0 2 since x\ > x 2 . Acting on the squared amplitudes yields: 

UYU 2 WX 2 + u\w 2 Wf 2 = - 0?0202 n 0i Ut >V]2 R , (2.23a) 

U 1 U 2 WX 2 + U 1 U 2 WX 2 = ~ ©^©IfW™ . (2.23b) 

The sum of these terms then gives: 

X^2W& = - ©l©2©2 n {^12 + 0?"^) 

X 

= - 0 p 0 2 0 2 n (-0? WfW$ + 0° ut Wi 2 r ) , (2.24) 


where we used the squared amplitudes from eqs. (2.20a) and (2.20b). Substituting into the 
expression of £ 2 (p) we obtain: 

£ 2 (p) = [ d^i©?©) 11 Wf x d<h 2 0£@!? W 2 r - I dn i2 @° ut @ 2 n W RR , (2.25) 

Jx 1>X2 j X\>X2 

where we introduced the shorthand notation dlli 2 ... m = n)=i d*!?*©?. In the hrst integral in 
the right-hand-side of eq. (2.25), the integrand is symmetric under the exchange k± -H- A; 2 , 
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which means that we can relax the condition x\ > X 2 and divide the integral by a factor 
2!. Hence this integral factors out into the product of two separate identical contributions 
from gluons k\ and k 2 , which both have exactly the same form as the one-loop result (eq. 
(2.15)). Thus we find for this term: 

Sf(p) = I d$i©?©i n >Vf x d$ 2 ©2©2 n W 2 R 

J X\>X2 

=% (-/ d^0?0i n Wf) 

=i(-C F a s L 2 ) 2 = i(sP) 2 , (2.26) 

which is just the expansion of the Sudakov E p (p) at second order. Hence: 


£ 2 (p) = £ p (p) + ^ G Go), (2.27a) 

^2 G (P) = ~ [ dn 12 0' i ,ut 0j, n W“ l . (2.27b) 

J X 1>X2 

The latter expression is the pure non-global contribution at this order. It is given by: 

S 2 G (p) = -^C F C A d 2 f -—dr ]1 dri2^-^ L e(x2-p)@(-Vi)Q(V2)A 1 q 2 g, (2.28) 

where we have 02 = 0(x 2 e _,?2 — p) k, 0(x 2 — p), since no collinear (double logarithms) are 
present for the pure non-global contribution. Hence the x integration easily factors from 
the rapidity integration and we just set the lower limit on x 2 to p (to single logarithmic 
accuracy). Performing the trivial integration over Xi we obtain the result L 2 / 2!. We note 
that at n th order we have: 

r 1 dai f xi dx 2 f X2 dx 3 

J P xi J p x 2 J p x 3 

Using the result of integration over cf >2 and 772 from eq. (A.6a) of appendix A (with { j, m} —>• 
{1, 2}), we obtain: 

= - ^CfCa|( f_j m 41n T -l 57 

= - = -|(c f CaC2 , (2.30) 

where L = a s L and £ is the Riemann-Zeta function. This is exactly the result obtained by 
Dasgupta and Salam [24] for NGLs at two-loops. To the best of our knowledge the analytical 
calculation of NGLs at finite N c beyond this order has not been performed before, and it 
is this very task that we do in the next section for the first time in the literature. 


rXn-l 

J p 


dx r , 


U 


71! 


(2.29) 


3 Non-global logarithms beyond leading order 
3.1 Three-loops calculation 

Having reproduced the well-known result for NGLs at leading order (two-loops), we proceed 
to compute NGLs at finite N c at next-to-leading order, namely triple gluons emission. As 
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usual we begin by the simplification of the measurement operator which will help us identify 
both angular and real-virtual configurations giving rise to large logarithms: 


U 3 = U.U2US = (l - 0?©i n 0f) (l - 0 P ©if0 R ) (l - 03©3 n 0f) 

= W 3 - Oi&2@3@3 n @3 (&2 + ©2 Ut ©2) (©1 + ©i"*©?) , (3.1) 


where U 3 is the collection of terms which when acting on the squared amplitudes W’^ 3 , 
with X summed over, yields a zero. The action of the measurement operator on the various 
squared-amplitudes summed over X gives: 

X] WaWfaa = + ©£ ut + ©? ut Wf 2 v 3 R + 0? ut ©2 Ut Wf 2 R 3 R ) . (3.2) 

X 

As stated in the introduction, the explicit expressions for the various squared amplitudes 
above (together with those at higher loops) will be presented in our forthcoming work [39]. 
Here we restrict ourselves to showing the simplification of the above squared-amplitudes in 
terms of the antenna functions defined previously in eq. (2.2). We have: 

YfizWva = -©I02©303 n x (©i n ©2 n WfW R >V R - ©i n 0^ ut w r W 23 r - 

X 

- ©/ ut 0/ n w 2 R W R 3 R - ©° ut © 2 n nf W RR + 

+©r W r J 3 r + 0 ? ut © 2 Ut W r 2 rr ) . ( 3 . 3 ) 


In eq. (3.3) W R and are defined at previous orders (see eqs. (2.12) and (2.21)), and 
W R 2 V 3 R and W R 2 R 3 R are the new irreducible terms of the squared amplitudes at this loop order 
proportional to the colour factor CfC^. 

Thus the hemisphere mass distribution at O(a^), S 3 (p), may be written as a sum of 
three contributions: £ 3 (p) = T, 3 (p) + £^(p) + T, 3 (p). The first contribution is: 

E$(p) = - [ d$i© 3 ©i n Wf x d$ 2 ©^0? W R x d$ 3 0^0“ W 3 r 

J Xl>X2>X3 

J d^eferwfj 3 = ^ (sf ) 3 = z *( P ), (3.4) 

which is simply the expansion of the Sudakov X p (p) at this order. The factor 1/3! accounts 
for the fact that the integrand in eq. (3.4) is completely symmetric under the exchange 
of gluons, which means that the condition x\ > x 2 > X 3 can be relaxed and the result 
multiplied by 1/3!. The integral is then factored out into the product of three identical 
integrals, each of them resembling the one-loop result eq. (2.15). 

The second contribution to £ 3 (p) is: 


Ef {p)=l dni 23 0 \ n w R 0 !) ut 0 i/w 2 3 ~+ / dn 12 30?w 2 R 0r©3 n w/3 ri + 

J Xl>X2>X3 j Xl>X2>X3 


tRR 


+ [ dn 1 2 . 30 :/w : / 0 r© 2 n >v; 2 \ 

J X\>X2>X3 


(3.5) 
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By swapping k\ -H> A: 2 in the second integral of eq. (3.5), and performing the successive 
permutations: k\ /c 2 then k\ k 3 in the third integral of the same equation, X^(p) 
becomes: 


sf (p) = [ dn 12 30 1 1 n wf©r© 3 n w^+/ dn 123 ©i n wf©r ©“ 

J X\ 


Xi>X2>X3 Jx 2>Xi>X3 

RR 


+ f dn 123 ©i n Wf 0 ^ ut 0 ^ W 23 

J X2>XZ>X\ 


(3.6) 


The three integrands in eq. (3.6) are identical except for the region of integration over 
transverse momenta fractions. Thus we unify them into a single integral with the region of 
integration expressed by: 

@(aq - x 2 )0(x 2 - x 3 ) + 0(x 2 — aq)0(a;i — x 3 ) + 0(x 2 - x 3 )0(x 3 — aq) = 0(x 2 - x 3 ). 

(3.7) 


Hence we write: 


S?(p) = ( - / dd>,0( , 0 i 1 n W| i 

=E?(p)xE» rG (p)- 


x - 


dn 23 ©° ut 0 |, n >v 


23 


' X2>X3 


(3.8) 


Thus X^(/?) factors out into a product of the one-loop primary cross-section and the two- 
loop NGLs cross-section (eqs. (2.15) and (2.27b) respectively). This result is expected 
from the expansion of the Sudakov form factor at one-loop times the leading NGLs. Said 
differently, S^(p) is just an “interference” term related to previous orders. 

The remaining term which is the pure irreducible NGLs contribution at this 

order, £ 3 G (p), is proportional to CfC^ and given by: 


S£(p) = ^ G (p) =~ [ dn 123 0° ut 0if (w RV 3 R + © 2 ut w RRR ) • (3.9) 

Jx 1>X2>X3 V ' 


From the above equation one sees that the new irreducible NGLs contribution at three-loop 
order is generated by two mechanisms: 

(a) the energy-ordered real gluons k\ and fc 2 outside Hr coherently emit the softest gluon 

k 3 into 'Hr — the term , 

(b) the hardest real gluon k\ outside Hr emits the softest gluon k 3 inside, while &: 2 is virtual 
(inside or outside Hr) — the term W R .^ R . 

In both cases NGLs result from the miscancellation with the corresponding squared ampli¬ 
tude for k 3 virtual, i.e., the miscancellation between and on the one hand, and 

_RVV ——RVR 

between W 123 and W 123 on the other hand. Both of these contributions are not related to 
previous orders. It seems, at first inspection, that the second mechanism mentioned above 
(particularly the case where gluon fc 2 is inside Hr, as clearly shown in eq. (3.10) below) 
is in contradiction with the common picture about the origin of NGLs. The latter picture 
dictates that, to all-orders, NGLs are entirely generated from a soft emission into Hr that is 
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coherently radiated by arbitrary ensembles of soft, but harder, large-angle energy-ordered 
gluons outside Hr [24, 25]. Nonetheless, and though NGLs contribution from the said 
mechanism comes from both gluons and k% inside Hr, gluon &2 is actually virtual. We 
shall see later that this mechanism persists at higher loops too. Hence whenever a contri¬ 
bution to NGLs comes from configurations whereby gluons other than the softest are inside 
Hr, then these “other” gluons must be virtual . 6 7 

In fact the contributions of the two terms W'\ 2 ;i and W 323 in eq. (3.9) are separately 
divergent but their sum is finite. The integral (3.9) can be expressed as a sum of two finite 
terms in the following way: 


V NG 

^3 


(P)=- 


L 


dn 12 30r©r 


Xl>X2>X3 


01 


rrrRVR . 

Wi 23 +0 


out 

2 


tttRVR rrrRRR 

W 123 +W 


123 


(3.10) 


Substituting the explicit expressions of the irreducible terms W 12 3 and W 12 3 in terms of 
the antenna functions yields: 


iNG 


(P) 


=^ CfC a J dr) 1 81n 2 (l - e 2l?1 )- 

- ^C F C| j drji 2 (A 2 q l(r)i) + Af q {r) 1 ) - 2C 2 ) , 


(3.11) 


where we performed the trivial integration over transverse momenta fractions to obtain 
L 3 / 3! and used the results of rapidity and azimuthal integrations shown in appendix A. The 
terms A 23 (r/i) and Af^r) 1 ) are given in eqs. (A.7b) and (A.7c) with {i,j, m} —>• {1,2,3}. 
The integration over rji in the first line of eq. (3.11) yields the result 8(^3 and that in the 
second line gives 4 ^ 3 . Thus the pure non-global contribution at this order reads: 

S 3 G (p) =^fC F CiC 3 - (3.12) 

Hence, up to this order we have: 


£ NG (p) =1 - ^C f CaC 2 + ^C F CiC 3 + 0(a 4 s ). (3.13) 

It is intriguing to note that the coefficients of NGLs at finite N c for two and three-loops 
(i.e., £2 and £ 3 ) are identical to those found at large N c [26, 30].' This is due to the fact 
that the combination of real/virtual squared amplitudes strangely produces identical inte¬ 
grands in the expressions of the NGLs contributions and (eqs. (2.27b) and (3.10) 
respectively). It would have been tremendously easy to resum NGLs for the hemisphere 
mass distribution to all-orders if the pattern in (3.13) persisted at higher loops. Although, 
we shall encounter Zeta functions at higher loop orders, the pattern itself unfortunately 
breaks down starting at four-loops, as we shall see in the next subsection. 

6 This observation was also made in ref. [30]. 

7 The appearance of £2 and £3 at two- and three-loop orders for NGLs is intriguing too. In appendix B 
we present some useful observations regarding possible relations between the two (NGLs and Zeta function) 
seemingly distinct quantities. 
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3.2 Four-loops calculation 

For four gluons emission, the computation of NGLs for the hemisphere mass distribution 
proceeds in an analogous manner to that of two and three gluons emission. At this (four- 
loops) order the measurement operator reads: 

U 4 = U4U2U3U4 = (l - Q?©”©*) (l - ©£©“ 0 *) (l - ©^©^©f) (l - ©^©j) 

= W 4 - 0?0^0^02©i n B4 (©3 + 03 ut @s) (©2 + ©2 Ut ©2) (©1 + ©r©f) , ( 3 . 14 ) 

where IA 4 is the sum of all terms which when operate on the squared amplitudes W’j ^ 34 and 
X is summed over give zero. Acting by the measurement operator on the various squared 
amplitudes and summing over configurations we obtain: 

X> 4 wf 23 4 = -©?0£0£©£©i n (WEST + ©l Ut Wi r 2 v 34 r + ©r W^ R + 

X 

+©3 ut +©r©r wi RR r+©r©^ wi v 2 r 3 rr + 

+ ©out©o ut yyRVRR + 0 o u t 0 out 0 out yyRRRR) (3.1.5) 

From eqs. (3.2) and (3.15), it should be clear how the result of the action of the measurement 
operator on the squared amplitudes at m th loop order would look like: 

• the softest gluon is always inside T~Lr , 

• each real gluon ki is associated with a step function 0 ° ut , 

• virtual gluons are associated with neither 0 in nor © out (i.e., they can either be in or 
out of T~Lr). 

The hemisphere mass distribution at fourth order may then be cast in the form (2.5). 
Substituting the various matrix-elements squared we can split the hemisphere mass distri¬ 
bution, £4 (p), into five parts: £4 = + £^ + £ 4 " + £^ + , with: 


Sf = 


dn 1234 e 1 ! 11 ©^©^©^ 1 , 


' Xi>X2>X3>X4 


Ef = - 


dni^fer ©?©?©! 11 w r 2 r w 3 r w 4 r + 2^3 + 204 + 


I Xl>X2>X3>X4 


£? = 


Sf = 


I Xi>X2>X3>X4 


' XI >X2 >X 3 >14 


+ [lf>3 and 2 «->• 4] + [2 -H- 1 then 1 o 4] + [2 1 then 1 3]^ 

dn 1234 {©? ut e“ 0 “ W 4 r (W R 2 V 3 R + ©r W RR 3 R ) + 3f> 4 + 

+ [3 2 then 2 H 4] + [1 f> 2 then 1 <4 3 then 1 <->• 

dni 234 ( 0 r 0 ? 0 r©i n W RR W R 4 R + 2h 3 + 2« 3 then 3 o 4 ) , 


Sf = - 


dn 1234 ©re^x 


' X\>X2>X3>X4 


1 pvout /-\out xaiRRVR x.out^out , a; RRRR\ /q i/?\ 

X (^1234 +©3 <^1234 + ©2 >^1234 + ©2 ©3 >^1234 J, (3.16) 
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where terms in the last line are the four-loop irreducible components of the squared am¬ 
plitudes for the corresponding gluon configurations. The parts Yf, E^, Yf, and Yf com¬ 
pletely reduce to integrals we calculated at previous orders, while the remaining Yf part 
is the new NGLs contribution. Let us evaluate each of these integrals separately starting 
with the reducible parts. 


3.2.1 Reducible parts 

For the first part E 4 we can, as usual, relax the condition x\ > X 2 > £3 > x 4 and multiply 
the result by a factor of 1/4!. This part then factors out into the product of four identical 
integrals of the form we met at 0(a s ) (eqs. (2.15) and (2.16)), thus we obtain for this term: 

Sf(p)=^(-/ d ^ i 0 ? 0 i n Wi R ) = i(E p ) 4 = E p (p), (3.17) 

which is just the expansion of the Sudakov at the fourth order. 

The second part Yf is carried out in a fashion analogous to that of E^ in eq. (3.5). The 
five integrands are transformed into the first integrand, , with the appropriate 

changes in the integration limits. Thus we have six integrals having identical integrands 
but different regions of integrations. Writing these integration regions as step functions 
and simplifying we obtain 0(xi — £ 2 ) ©(2:3 — X4 ). This means that Y f factors out into the 
product of two integrals, one over k\ and A : 2 and the other over k 3 and Aq, as follows: 

sf (p) = ( I d<f> 3 @ p 3 ofw$ x x (- [ dni20? ut 02 n >Vi2 

VI 3>I4 / V J X\>X2 

=|(Ef ) 2 x S* G , (3.18) 

where we have used eqs. (2.26) and (2.27b) to arrive at the second line of the above 
equation. Eq. (3.18) is in fact the interference of the expansion of the Sudakov form factor 
E p with NGLs at two-loops E^g. 

Performing the integrations in the third part Yf along the same lines outlined above 
for Yf (i.e., by making appropriate changes of variables) we obtain: 


YX(P) = [- ^ 40 ^ 0^4 * ~ 


' X1>X2>X3 


dn :L 230r t 03 n 


tttRVR 

W 123 +e 


,out 


w 


;RRR 

123 


=e Up) x E^G (P ) 


(3.19) 


where we used eqs. (2.15) and (3.9). This result is the interference between the expansion 
of the Sudakov and Eg G . 

The part Yf can be written as follows: 

Yf(p)=l[ dn 12 0? ut 02 n Wi^ x f dn 3 40 r 04 n W*r 

^ Jx 1 >X 2 JX3>X4 

=i(ENG) 2 , ( 3 . 20 ) 

which indicates a possible pattern of exponentiation of NGLs since this term resembles the 
structure of the expansion of exp{E^ G | a t this (fourth) order. 
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3.2.2 Irreducible part 

The irreducible part at four-loops, Tif(p), is given by: 

Ef ( p )=-[ dn 1234 0r©i n x 

J a;i>a;2>3:3>^4 

x (w rv 34 r +© r w rvrr +© r w rr 4 r +© r©r w rrr 4 r ) . ( 3 . 21 ) 

The first irreducible squared-amplitudes and W R 234 R are proportional to CfC^, 

whereas the other two amplitudes, W R 234 and W] 1 ™, contain both CfC^ and CpC^ 
terms. The phase space integration of all terms in eq. (3.21) breaks the simple pattern 
observed in eq. (3.13), with the last two terms (Wi 234 R and W 1234 ) even breaking the 
colour pattern. Therefore the loss of the pattern of NGLs is in fact a manifestation of the 
break in the structure of the eikonal amplitudes. This might be related to the failure of the 
“probabilistic scheme” discussed by Dokshitzer et al. in ref. [32], where such (irreducible) 
contributions to the eikonal squared-amplitude were dubbed “monster” terms. They were 
traced back to be originating from the “colour polarisability” of jets [32]. 8 Moreover, we 
note here, as can be seen from the form of , that contributions to NGLs at this order 
are generated when the softest gluon £: 4 is emitted inside the measured hemisphere HLr, 
whilst the hardest gluon k\ is always outside. The other two gluons, & 2 and fc 3 , may be 
emitted inside HLr provided they are virtual (in accordance with the observation made in 
the previous subsection). 

The contribution Ylf is not related to the expansion of the Sudakov nor to NGLs at 
previous orders, and represents the new non-global contribution at four-loops. Together 
with they form the total non-global contribution, £4 , to the hemisphere mass distri¬ 
bution at this order. To evaluate the part we first integrate over transverse momenta 
fractions, which as usual yields L 4 /4!, and then perform the azimuthal and rapidity in¬ 
tegrations according to the angular configurations indicated in eq. (3.21). The azimuthal 
integrations are not as straightforward as at three-loops and we find the method used in ref. 
[30] of contour integration very useful in reducing the number of integrals to be performed. 
The reduced integrals are then carried out analytically whenever possible, otherwise numer¬ 
ically. Since all analytical integrations yield results that are explicitly proportional to £ 4 , 
the resultant values from numerical integrations were interpreted in terms of £ 4 . In addition 
to this semi-numerical approach we also verify our results by numerically integrating each 
of the finite terms in the following form of T,f(p): 


sf (p) = - diwereirx 

J X\>X2>X3>X4, 

x {ei?ei i "W“ v ™+ 

+ ejfer (W™ + w™) + ere? (W?™ + W“ 3 V 4 ") + 

+er'er (w™+w“™+w™+w““)}, (3.22) 

8 More details are to be found in our forthcoming paper [39]. 
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over the full (7-dimensional) phase space using the multi-dimensional numerical-integration 
library Cuba [45]. The final result reads: 


Sf (p) = - 


L 4 / 25 


4! V 8 


-c F ci + c|cl) c 4 , 


which may also be rewritten in the following two alternative forms: 


sf(p) = -^c F cic 4 


29 

— + 


Cf 

Ca 


Z 4 

" 4! 


25 o 2 9 

— C F C| C 4 + jr (C f Ca C 2 ) 


(3.23) 


(3.24a) 

(3.24b) 


The expression (3.24a) explicitly shows the finite- N c correction to the large- N c result, while 
(3.24b) emphasises the pattern C f Ca Cn+i seen at two and three-loops. It also reveals that 
even though is a new irreducible contribution at four-loops, it still contains factors 
related to lower-order NGL contributions (the term (C f CaC 2 ) 2 )- Observe that the size of 
the finite- N c correction in (3.24a) is about ~ 1.5% that of the large- N c result. This is in 
agreement with the conclusion arrived at in [46] for the impact of finite- N c corrections at 
all-orders for e + e _ processes. The total non-global contribution at this order is then given 
by: 


L 4 / 25 


£4 °(p) = S f + Ef = ~ ^ C F CiC 4 - ^ C 4 Ci Cl 


13 


(3.25) 


and thus the hemisphere mass distribution up to this order is expressed as: 


£(p) =£ F (p)xS NG (p), 


(3.26) 


iNG 


L 


L 3 


L 4 


(, P) =1 - ^c f CaC 2 + —c F cic 3 - — 


3! 


4! 


fc F c|c 4 -f c|cicl +o(«: 


In the next subsection we discuss the five-loops case and the possibility of resummation of 
NGLs. 


3.3 Five-loops and beyond 

Following the same steps as before we write the measurement operator at five-loops as 
follows: 

W 5 = (l - 0?©i n 0f) (l - 0£0!f0£) (l - ©W©!) (l - ©W©") (l - 0£©?©f) 

=U 5 - 0f0^0g020g©^©^ (el + ©^© 4 ) (©3 + ©3 Ut 0 3 ) (©2 + © 2 Ut © 2 ) X 

x (©^ + 0° ut @^) , (3.27) 

where, as usual, is the sum of all terms that yield vanishing contributions to the hemi¬ 
sphere mass distribution. Acting by the measurement operator on the various squared 
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amplitudes and summing over configurations we obtain: 


E^Wf 23 45 = -n 0 i 0 5 


X 


(W| 


i=l 


iVVWR i /'AOUt-iAiRVWR i /'AOUt-i aiVRWR , pvORtlA jVVRVR , /-\OUt-i aiVWRR , 
x 12345 + U 1 ' / ' / 12345 + U 2 |/ ' / 12345 + u 3 Kl/ 12345 + u 4 |/ ' / 12345 + 

+ 01 Ut 02 >Ut W[ 2345 R + 0° ut 0° ut Wf 2 ™ + 0 OUt 0 OUt W RVVRR + ©OUtQOUtyyVRRVR. 
+ 00^00^-WVRVRR + 0 out QOUtyy WRRR - — - — - — ~ 


’ IZO^O 1 1 4 ' r 1 ^ 

dVVRRR I rxoutr^out/^vout^ AlRRRVR 

■ ~z ' izc^to 1 ' ^12345 “T w l w 2 w 3 KV/ 12345 

+ ©^©^©OUly^RRVRR + QOUt QOUtQOUty^RVRRR + 0 out 0 out 0 out 
I r\ out r\ out r\ out out ;RRRRR \ 


The hemisphere mass distribution at five-loops is then given in eq. (2.5) with m = 5. 
Substituting the various matrix-elements squared and following the same procedure outlined 
at four-loops we again obtain two types of contributions; reducible, ££, and irreducible, £g r . 
The former contains all the interference terms between the Sudakov factor £ p and NGLs at 
previous orders as well as interference terms between two and three-loops NGLs. Explicitly 
written it reads: 


-£(**)■ + M 


(Ef) 3 x E? g + 


2 ! 




P\ 2 


x + £ p x E. 


NG + SfxSf 


(3.29) 


Note that the penultimate term in the above equation contains the contribution E p x 
(£^ g ) 2 /2! (eqs. (3.25) and (3.20)). 

The irreducible contribution E~ rr is expressed as: 


virr 

^5 


dn 1234 5 ©re in ' 


' X\>X2>X$>X4 : >X5 


tttRVVVR , 

Wl 2345 +© 


° ut v\7 KKV 
2 * / * / 12345 


+ 0|| ut 


. . .RVRVR , ^, 
^12345 +©: 


° ut v\7 H VK 
4 Kv 12345 


+ 0^ Ut 


0 ° Ut 


■» . ,-RRRVR 
^12345 + 


, rNOUtpvOUt^ a i rNOUtrNOUtA a ,RvRRr 

-hU 2 U 4 KV 12 345 -h ^3 KV 123 45 


+ © 2 ut 


© out /^out u ,RRRRR \ 

3 U 4 KV 12345 J , 


(3.30) 


which is neither related to the Sudakov factor nor to NGLs at previous orders, and which 
contains both CpC^ and CfC^ terms. The irreducible squared-amplitudes that contribute 
to Sg rr can be classified into two types: 


, . , . , „ „ 4 tttRVRRR t-ttRVRVR ^tttRVVRR . prrRVWR 

proportional only to CpC^: W 12345 , W 12M5 , >V 12345 , and W 12345 , 


. , , -a ^ ^ 4 t-ttRRRRR tttRRRVR 

• containing both and CpO^: n / 12345 , W 12345 , 


w 


12345 


and W 


RRVVR 
12345 > 


which 


all contain “monster” terms. 


The calculation of Eg r turns out to be trickier and more involved than anticipated. In 
particular we have not yet been able to simplify (like we did at four-loops) the monster 
parts of the irreducible amplitudes of the second type above, to forms that can readily be 
integrated. Other than the monster parts, all remaining terms (either of the first or second 
type above) are in fact integrable. 
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Till the full expression of the irreducible contribution is simplified and integrated, we 
express the result of £ij rr in the following form (based on the pattern seen at four-loops 
(3.24b) and the pieces found in the integrable amplitudes of eq. (3.30)): 


V'irr 

^5 


^CpClCs 


ol T /5 


/Cf _ 1 

Ua 2 



C F ClC 5 + a/3C|CiC2C3 , 


(3.31a) 

(3.31b) 


where a = C5/C2C3 — 0.5244 and the constant coefficients a and /3 are yet to be determined. 
We discuss the possible values of these constants when we compare our results with those 
at large N c in the next section. The form (3.31a) explicitly shows the finite-A r c correction. 
The total NGLs contribution at five-loops then reads: 


V NG 

^5 


= Y^ G X s* G + Si; rr 

= -^T c FGiC2C3 + si i r 


5! 



C F CaCs — (10 — a/3) C 2 f C 3 a ( 2 C 3 


(3.32) 


The results we obtained up to five-loops, particularly eq. (3.29), in fact suggest a 
possible resummation of NGLs into an exponential function of the form: 


£ NG (p) = exp <| - |^C f CaC 2 + ^C F CiC 3 - ^ C F Ci C 4 


29 

¥ + 


Cf 

C a 


+ ^c F clc 5 


ol + (3 


Cf 

Ca 


+ 0{al) 


+ 


(3.33) 


Eq. (3.33) may actually be rewritten in a form analogous to that found in ref. [47] (eqs. 
(5.10) and (5.11)) for clustering logarithms. To this end we write: 


S NG (p) = exp 


— V —5 

C A ^ n\ Sn 

n>2 


(-C A I ) 1 


(3.34) 


where 


$2 — ( 2 , 


S3 — C3 , S 4 





55 


Cs 




(3.35) 


The above-mentioned similarity between clustering logarithms and NGLs emphasises the 
common (non-global) origin of the two types of logarithms. Moreover, following the pattern 
in (3.31b), eq. (3.33) may also be recast into the form: 


X NG (p) = exp \ - ^C f CaC 2 + ^C F CiC 3 - ^ 


fc F ciC4 + ^c|c| cl 


+ 


+ 


5! 


a — 


P 


C F ClC5 + a/3C2c|C2C 3 


+ 0 {a\ 


(3.36) 
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The expansion of the above exponential exactly reproduces the terms we have calculated 
up to five-loops including all interference terms in the distribution. At each higher order 
one simply adds a new irreducible NGLs term in the exponent. 

In fact, if the pattern deduced in eqs. (3.24a) and (3.24b) persists at higher-loop orders, 
then one can put forth the following ansatz for the general form of the n th order contribution 
to the exponent of the resummed NGLs factor: 


T n 

(—l)- 1 GpCr'Cr. 

n\ 


(-i) 


n—1 


L n 


n\ 


Ln/2j 

7 n + 0 " k 

k =2 

Ln/2j 


C F 

Ca 


jfc-1 


2 k ~i 


7n C F Ca _ 1 Cn + ^ ^ Cl ~ k CfcC n-k 


k =2 


(3.37a) 

(3.37b) 


where \n\ represents the floor function of n, and 7 „, a n , 7 n , and are constant coefficients 
to be determined from integrations. The two formulae presented above for the ansatz are 
equivalent up to the constant coefficients. The first form stresses the finite-iV c correction 
while the second preserves the pattern seen at two-, three- and four-loops (eq. (3.24b)). 
The above formulae may only be verified once higher-loop orders are carried out explicitly. 
We hope to perform such calculations in the near future. 

We note that in the exponent of (3.33) (and (3.36)) the series in L has alternating signs 
at each escalating order. To assess the relative size of the three and four-loops corrections 
to the leading two-loops result, we plot in figure 2 the ratio E NG /exp(£^ G ) for various 
truncations of the series in the exponent in eq. (3.33). The leading NGLs coefficient seems 
to dominate for only relatively small values of L (L < 0.15). For larger values the series 
seems to depart from the leading term in an alternating way (towards larger (smaller) values 
for odd (even) loop orders). These significant variations mean that the terms computed thus 
far are insufficient to capture the full behaviour of the all-orders resummed distribution. 
We expect, however, that adding few more terms in the exponent may lead to a convergent 
and more stable behaviour, since one could argue that higher-order terms are suppressed 
by L n /n\, while saturates at 1 as n becomes larger. 

To cross-check our results, eqs. (3.33) and (3.36), we compare them, in the next section, 
to previous calculations at large N c both at fixed order and to all-orders. 


4 Comparison with large- N c results 

4.1 Comparison with analytical results at large N c 

Having calculated the coefficients of NGLs at finite N c fully up to four-loops and partially 
at five-loops, we can now compare our findings to those of Schwartz and Zhu [30] obtained 
through the analytical solution to the BMS equation [26] in the large-iV c limit. To go 
from the finite-iV c case to the large-A r c approximation we simply invoke the replacement 
C F —>• Ca/2 = N c /2 (where Ca = N c ). This is equivalent to expanding C F to first order in 
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Figure 2. Plot of the ratio S NG (p)/exp(E^ G (p)) in terms of the logarithm L = a s / 7 rln(l /p). 


colour: 


C F 


N 2 C ~ 1 
2 N c 



(4.1) 


Then the expansion of our result (3.36) to leading order in colour, i.e., at large N c , up to 
five-loops is: 


S NG =1 - —(N C L ) 2 + —(N c Lf + 

24 v c ) t 12 v c ) v- 34560 

7r2 C 3 . ( 2 a — /3)Cs + 


7r 


+ - 


288 


+ 


480 


(N C L) 4 + 

I {N c Lf + O ((. N C L) 6 ) , (4.2) 


where we have written the explicit values of £2 = 7 t 2 /6 and £4 = 7t 4 /90. The full result 
reported by Schwartz and Zhu (SZ) at large N Cl Sg G ) is [30]: 


V NG 

^SZ 


=1 - —L 2 + —L 3 + 
24 12 


7T 


34 560 


L 4 + 


V 360 480 


L 5 + 0(a 6 s ), 


(4.3) 


where L is simply N C L. The two results are thus identical up to four-loops. Recalling that 
they were arrived at using different approaches, their equality provides a solid cross-check of 
the correctness of the computed NGLs coefficients (at least up to four-loops). Even though 
we are unable to fully compare our result to that of Schwartz and Zhu at five-loops, due to 
the missing values of a and /3, it is ironic to note that the pattern spotted at four-loops, 
eq. (3.24b), seems to hold true at five-loops. Whilst the term £5 is apparent in (4.3), the 
product (( 2 C 3 i s disguising in the factor 7 t 2 <^3/360. A quick comparison reveals the values: 


17 1 

a — — + - 
2 a 



(4.4) 
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Given the above values we expect the finite-A r c result of NGLs up to five-loops to be 
expressed as: 


S NG (p) = exp \ - ^C f CaC 2 + ^CpCiCs - ^ 


+ 


L 5 

5!~ 


y C F Cl C5 + 2 C^C| C2C3 


f c F c|C4 + ^c2cic 2 2 


+ 0 (a«) 


+ 


(4.5) 


We are not, however, claiming to have fully accounted for NGLs at this order since we have 
not explicitly calculated the coefficients of NGLs with colour factors CpC^ and C F C^. We 
hope that further research on this can help verify the above equation (and eq. (3.37) in 
general) in which case it may actually be possible to find a key to the analytical resummation 
of NGLs both at large and finite N c . 

Our approach additionally has the benefit that it sheds light on the possibility of 
assessing the validity of the large-lV c approximation, by means of judging the impact of 
neglected finite- N c corrections. An important note in this regard is that at two- and three- 
loops, as can be seen by comparing eqs. (3.26) and (4.2) at 0(L 2 ) and 0(Z 3 ), there are 
no hidden terms buried by the large- N c approximation, and the finite- N c result can simply 
be obtained from the solution of the BMS equation by just restoring the full colour factors 
through: N 2 —>• 2C f Ca (at two-loops) and IV 3 —>• 2C F C^ (at three-loops). At four-loops 
this is not true and in fact there is a hidden correction that is given plainly in (3.24a). 
We regard this as the first-order proper fmite-lV c correction which introduces new terms 
that are entirely absent at large N c . The second-order proper finite- N c correction occurs 
at five-loops and is shown in (3.31a). 


4.2 Comparison with all-orders numerical results 

In order to verify our resummed formula (3.33), and even (4.5) which includes the five- 
loops term, it is instructive to compare it to the all-orders numerical solution of the finite- 
N c Weigert equation [28]. We have not, unfortunately, been able to obtain the output of 
the MC program, written by Hatta and Ueda [27], for the hemisphere mass distribution. 9 
We thus postpone this discussion till the said numerical distribution becomes available. 
Furthermore, to assess the importance of the missing higher-loop terms in (3.33) we compare 
it to either the results obtained by the numerical solution of the BMS equation [26] or to 
the output of the numerical MC program of Dasgupta and Salarn (DS) [24], The latter two 
numerical solutions are in fact identical within a percent accuracy [26, 30], and we thus 
restrict ourselves to the DS MC program. 

Let us introduce the standard evolution parameter t [24], which accounts for the running 
of the coupling: 



1 

4vr/3 0 


In 


1 

1 ~ 2/?o u s L 


9 The hemisphere mass distribution has not yet been coded into the MC program [46]. 


(4.6) 
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Figure 3. Plot of the NGLs function E NG (p) at large (left) and finite (right) N c . 


where f3o is the one-loop coefficient of the QCD /3 function. At fixed order one has t = 
a s L/2ir = L/2. Hence substituting L by 2t into eq. (3.33) up to four-loops we find: 


£ NG (f) = exp ( -C F C A y t 2 + ^C F C|C 3 1 3 - ^ ' 


25. 

8 " 


rd r^'6 I r^2 r\2 


f 4 + 0(f 5 ) 


(4.7) 


We compare the result (4.7) with the parametrisation for NGLs to all-orders obtained in 
ref. [24] by fitting to the output of the aforementioned DS MC program [24]: 


£i)S (*) =ex P 


(-C F C A 
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1 + (0.85C A i) 2 
1 + (0.86C A f ) L33 



(4.8) 


In figure 3 we plot our approximate resummed result (4.7) for various truncations along 
with the DS resummed factor (4.8) for a range of t € [0,0.5] both at large (left) and finite 
(right) N c . Recall that a value of t = 0.3 corresponds to a value of L = 19 and p ~ 10 ~ 8 
for a s ~ 0.1, which is sufficient for phenomenological purposes. Few points to note from 
the plots. Firstly, as expected for finite N c , all curves are shifted up due to the fact that 
one is using C F = 4/3 ~ 1.33 instead of C F = C A /2 = 3/2 = 1.5 (recall that C F is in the 
exponent). Secondly, it is striking to observe that the best approximation to the all-orders 
result for quite a large range of t is the leading two-loops result , 10 exp (S XG ), for both large 
and finite-Ay cases. This suggests that the alternating, positive and negative, higher-loop 
contributions to X NG somehow balance out. 

Moreover, the main feature of the plots and which has direct link to the purpose of 
this paper is actually seen at small values of t. We see that the interval of t over which the 
four-loops result and the all-orders resummed factor overlap is 0 < t < 0.12. This interval 
of overlapping is smaller for three-loops, 0 < t < 0.08, and even smallest for two-loops 
0 < t < 0.05. The latter feature may be seen more clearly in figures 4 and 5. One would 
therefore expect that adding more terms in the exponent of (3.33) leads to increasingly 
larger intervals of overlapping. 

10 This observation was made in refs. [29, 30] too. 
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Figure 4. Plot of the ratio S NG (p)/Epg 
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A similar observation was also made in ref. [29] for “filtering analysis” in the case of the 
filtering parameters ngit = 2 and r/fu t = 0.1, 0.3 (figure 18 of ref. [29]), as well as ngi t = 3 and 
^filt = 0.3 (figure 21). However the author of ref. [29], having plotted the expansion of the 
full filtered Higgs-jet mass distribution including both primary and non-global logarithms 
in the Cambridge-Aachen jet algorithm [48, 49], ascribed the convergence of the series, 
as one adds higher-loop terms, to the dominance of the primary series. The author then 
verified this explanation by plotting the same distribution for higher values of r/ gi t (figure 
19) where collinear logarithms are expected to be absent and NGLs become of the same 
order as primary logarithms. Two issues to point out regarding our work compared to that 
of ref. [29]: firstly, we are plotting purely the NGLs resummed exponential factor and 
hence the convergence seen in figures 4 and 5 has nothing to do with primary logarithms. 
Secondly, it is well known [43, 50-53] that employing the Cambridge-Aachen jet algorithm 
not only reduces the size of NGLs but also introduces clustering logarithms that are as 
important as NGLs. Thus plotting the full distribution, which includes primary, non-global 
and clustering logarithms, would not tell much about the convergence of the NGLs series. 

Moreover, the author of ref. [29] also plotted (figure 24) the pure NGLs resummed 
factor for the interjet energy flow distribution and concluded that, up to six-loops, the 
NGLs series seems to be divergent. Recalling that the coefficient of the two-loops NGLs 
depends on the rapidity gap At/ [25, 43], it is likely that higher-loop NGLs coefficients 
depend on Ar/ too. The divergence may thus be due to the presence of the A rj terms. 
For the hemisphere mass distribution that we have treated in this paper there is no such 
rapidity gap dependence. A proper answer, however, may only be given once the former — 
interjet energy flow — distribution is carefully considered, a task which we hope to perform 
in coming publications. 

Notice finally that the so far discussed NGLs behaviour is in contrast to that seen 
for clustering logarithms [47] where the whole structure of the all-orders result is mostly 
captured by the first few terms in the exponent. 
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Figure 5. Plot of the NGLs function E NG (f) at large N c for the range 0 < t < 0.14. 

5 Conclusions 

In this paper we have considered the calculation of the leading NGLs at finite N c up to 
five-loops for the hemisphere mass distribution in e + e~ —> di-jet events. We performed 
the calculation by means of integrating the squared amplitudes for the emission of energy- 
ordered soft gluons in the eikonal approximation, valid at single logarithmic accuracy, over 
a suitable phase space achieved through a measurement operator. The two and three-loops 
results were shown to be relatively straightforward to obtain, and were found to be directly 
related to the Riemann-Zeta function. We noticed that, up to this loop order, fmite-lV c 
corrections are absent. This was a direct consequence of the relatively simple structure 
of the eikonal amplitudes (up to this order) as well as the combination of real/virtual 
amplitudes induced by the phase space measurement operator. 

Within the same eikonal framework we computed the four-loops contribution to NGLs 
distribution. The latter turns out to be much harder than the previous two orders and 
the simple result in terms of a product of a single colour factor and a Zeta function breaks 
down. This failure originates from a break in the simple structure of the corresponding 
real-virtual eikonal amplitudes at four-loops order, a phenomenon that was noticed more 
than two decays ago [32]. Nevertheless, we were able to overcome this complexity, compute 
NGLs and even spot a new pattern for NGLs at and beyond this order. This pattern helped 
us to successfully write down the five-loops contribution to NGLs up to constant coefficients 
which we extracted from comparisons to previous large- N c results [30] . We hope to be able 
to fully compute these constants in the near future. The five-loops calculation reveals that 
the NGLs distribution seems to exhibit a pattern of exponentiation. To this end the said 
distribution was cast in an exponential form with full NGLs coefficients and colour factors 
up to four-loops in the exponent. 

Comparisons to large- N c results obtained by other authors [30] confirmed our findings, 
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at least in the latter limit. We then took the step forth and compared our exponential 
function to the all-orders numerical resummed result reported by Dasgupta and Salam 
[24]. To our surprise, the shape of the all-orders result was best represented by the two- 
loops approximation for a wide range of the evolution parameter t. In the region of small 
t, however, adding more terms in the exponent of our resummed result yielded better 
agreement, then the two-loops result, with the all-orders numerical result. This suggests 
that more higher-loop contributions are needed for our result to be of any phenomenological 
significance (i.e., till the agreement extends to values of t up to ~ 0.2 - 0.3). The task of 
computing these higher-loop terms might not be impossible after all given that we have 
developed, over the course of preparing this paper, the machinery for: computing eikonal 
amplitudes at finite N c to theoretically any loop order, reducing the dimension of the phase 
space over which to integrate, and spotting a pattern for NGLs at each order. 
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A Angular integrations 

In this section we present some definitions and azimuthal/rapidity integrations which have 
proven useful in our calculations. 

Following ref. [30] we define the following angular functions: 


(■ ij) =cosh(? 7 i - rjj ) - cos (fa - fa ), 

Uj) = _ M _. 

2 sinh rji sinh rjj 

The basic antennas are expressed as: 
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The ^-azimuthal averaging over the inverse of the angular function (ij) is given by: 


(A.2) 


l 


27r dfa 1 ,, , 

__ = cschte - v \ 


(A.3) 


In the case where a gluon k m is constrained within the measured hemisphere region and 
gluons ki and kj are constrained outside, the k m angular integration over w™ yields [30]: 
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with r]i < 0 and r/j < 0. Furthermore, the azimuthal average over the angle 0j of the emitter 
ki yields [30]: 
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(A.5) 


With the same conditions (rji < 0, r/j < 0 and rj m > 0) we can perform the azimuthal 
and rapidity integrations for the following antenna functions: 
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We can also perform further integrations over cj)j and r/j: 
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(A.6c) 
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with Li 2 the poly logarithm function of order 2. Similarly integrating over the angles of the 
softest particle k n in the antenna A™ n yields: 


A=wYf (ln(l + (im)) + ln(l + (jm)) - ln(l + (ij))) 
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At four-loops, the following azimuthal integrations are relevant: 
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Then we have: 
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B A note on NGLs-£ n relation 


As a byproduct, we notice from eqs. (2.30) and (3.11) that one may define the following 
possibly “new” logarithmic-integral representation for the Riemann-Zeta function: 
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where T(s) = (s — 1)! is the Gamma function and the variable s is greater than 1. The 
above formula is valid if s is an integer. In the case of non-integer real values one has to 
take the modulus of the right-hand-side in eq. (B.l). In terms of the polar variables (0,(f)) 
the Zeta function admits the integral formula: 
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where c = cos0. Notice that the form (B.2) seems to fail (in Mathematica 9) for s > 10. 
Moreover, it is interesting to note that the first non-divergent value of the Zeta function 
is for s = 2, and so is the first non-vanishing coefficient of NGLs. If we let S s denote the 
NGLs coefficient at the s th loop order then we can write S s as the Mellin transform of the 
function (e 71 — l) -1 [54]. That is: 

1 r + °° 77 s " 1 

S s = Cs = TTTT / -d?7) (B- 3 ) 

s r(s) Jo e 71 - 1 

which is at least true for two and three-loops NGLs coefficients. Recall that the Mellin 
transform techniques were employed in ref. [55] to compute the first resummed result for 
event-shape distributions. Whether there exists a more profound relation between NGLs 
and the Zeta function (and its related functions such as the polylogarithms) is a subject 
that requires further investigations. 
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